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TECHNICAL FIELD 

The present invention is directed to a method of processing a 
signal, for example an optical coherence tomography signal, in the time 
domain and more particularly to such a signal processing method using 
non-linear regression with a sinusoidal model. 

BACKGROUND OF THE INVENTION 

There are a number of known techniques for imaging tissue 
including ultrasound tomography and optical coherence tomography. In 
ultrasound tomography, a sound wave is launched into tissue where it is 
reflected off of interfaces, the sound energy returning to a detector. 
These interfaces are places where characteristics of the sound energy 
changes. The intensity of the returning signal is plotted as a function of 
time, where the time axis is transformed into a depth axis using the 
knowledge of the approximate speed of sound in the tissue to obtain an 
intensity profile as a function of depth. The plot of intensity as a 
function of depth along a particular lateral location is known as an A- 
scan. From a collection of A-scans, each obtained from neighboring 
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lateral locations, an image depicting intensity as a function of lateral 
position and depth can be obtained. Such an image is called a B-scan 
and is commonly used in obstetrics to analyze a growing fetus, in 
cardiology, in neonatal encephalogy and in certain aspects of internal 
medicine. 

In optical coherence tomography, an optical wave is launched 
into tissue where it is reflected off of interfaces, the optical energy 
returning to a detector. The interfaces are places where the 
characteristics of the light, e.g., the complex index of refraction, 
changes. The intensity of the returning signal is plotted as a function of 
time, where the time axis is transformed into a depth axis using the 
knowledge of the approximate index of refraction in the tissue to obtain 
an intensity profile as a function of depth. A plot of the intensity as a 
function of depth along a particular lateral location is known as an A- 
scan. The major differences between ultrasound tomography and optical 
coherence tomography relate to the differences between sound and light. 
Sound is slower than light. Thus, in ultrasound tomography, to develop 
the image, the time between the emission of sound from the source and 
the time when the returning sound energy is collected by the detector 
can be measured using a clock. Because light travels so quickly, an 
interference technique is used to measure time. Thus, optical coherence 
tomography is based upon an interferometer. From a collection of the 
A-scans obtained via optical coherence tomography, where each A-scan 
is obtained from neighboring lateral locations, an image can be 
generated depicting intensity as a function of lateral position and depth. 
Again, such an image is called a B-scan. With known ultrasound 
tomography systems, the scan can interrogate on the order of 10cm into 
tissue with a spatial resolution of approximately lOO^im. With prior 
optical coherence tomography systems, the scan can interrogate on the 
order of 4mm into the tissue with a spatial resolution of approximately 
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l^m. The depth and resolution of both imaging techniques are 
dependent on the source, detector, and signal processing employed. 

The optical coherence tomography signal is typically processed 
in the frequency domain. However, the frequency domain is limited by 
5 the inverse duality of spatial resolution and frequency resolution. 

Specifically, the finer the spatial detail obtained, the coarser the 
frequency detail obtained and vice versa. Thus, with processing in the 
frequency domain, one cannot obtain both the high spatial resolution 
needed to see small vessels and the small frequency shifts needed to 
10 quantify the low fluid velocities in such vessels. In order to overcome 

this problem, known optical coherence tomography method have 
averaged a large number of A-scans, which requires expensive 
technology. 

1 5 BRIEF SUMMARY OF THE INVENTION 

In accordance with the present invention, the disadvantages of 
prior signal processing techniques for processing signals such as an 
optical coherence tomography signal, have been overcome. With the 

20 method of the present invention, both high spatial resolution and high 

frequency resolution are obtained with less sophisticated instrumentation 
so as to provide a lower cost system and method. Although, the 
invention is described below for processing an optical coherence 
tomography (OCT) signal, the method of the present invention is not 

25 limited thereto. The method is applicable to the processing of other 

signals where the frequency of the signal is in a known range. 

More particularly, the method of the present invention includes 
digitizing the signal to be processed to provide digital data points. The 
digital data points, representing a portion of the signal, are processed 

30 using non-linear regression with a sinusoidal model optimized for a 

known frequency range. The parameters or coefficients of the sine wave 
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fit to the digital data points provide the information for which the signal 
is processed. Specifically, in one embodiment of the present invention, 
the parameters or coefficients are used to generate high resolution A- 
scan and B-scan images depicting fluid flow. In addition to information 
relating to the velocity of fluid flow, shear rate information is obtained 
as well. 

In one embodiment of the present invention, the digital data 
points that are processed represent a portion of the signal that is less than 
a full cycle of a wave of the signal. As a result, higher spatial resolution 
is obtained. In accordance with the processing method of the present 
invention, a single A-scan can provide good images. The quality of the 
images can be improved by using only a few overlapping A-scans as 
opposed to averaging of a large number of A-scans as in the prior art. 
Thus, with the method of the present invention, turbulent flow in small 
vessels can be imaged from data collected in a short period of time. 

These and other advantages and novel features of the present 
invention, as well as details of an illustrated embodiment thereof, will be 
more fully understood from the following description and drawings. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS 
OF THE DRAWING 

Fig. 1 is a block diagram illustrating an OCT system providing a 
signal for processing in accordance with the present invention; 

Fig. 2 is a graph illustrating digital data points of a noisy sine 

wave; 

Fig. 3 is a graph illustrating a sine wave fit to the digital data 
points of Fig. 2 in accordance with the present invention; 

Fig. 4 is a plot of velocity versus depth obtained with the 
sinusoidal model of the present invention and no filtering; 

Fig. 5 is a plot of velocity versus depth obtained with the 
sinusoidal model of the present invention with median filtering; 
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Fig. 6 a plot of velocity versus depth obtained with the sinusoidal 
model of the present invention with the removal of points that have 
failed to converge accurately; and 

Fig. 7 a plot of velocity versus depth obtained with a normalized 
sinusoidal model of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

The signal processing system of the present invention focuses on 
a small section of a complex signal, such as an OCT signal, and models 
this small section of the signal as a single sine wave for time domain 
processing. This is opposed to prior methods using frequency domain 
processing of OCT signals that sample multiple sine waves of multiple 
frequencies, amplitudes, and phases. The time domain processing 
method of the present invention is described in detail below for an OCT 
signal, but the method is not limited thereto. The signal processing 
method of the present invention can be used for signals, the frequency of 
which is known within a given range. The size of the frequency range, 
within which the frequency sought is known to be, varies depending on 
the signal to noise ratio such that the higher the signal to noise ratio, the 
greater the frequency range can be. 

An OCT system that generates the OCT signal to be processed as 
shown in Fig. 1 is based on a Michelson interferometer. The 
interferometer includes a sample arm 10 and a reference arm 12. Light 
in the sample arm is reflected by the interfaces in the tissue sample 14 
and interferes with the light in the reference arm 12. When the optical 
path lengths of the reference and sample arms are equal to within the 
coherence length of the source 16, an interference signal is produced. 
The coherence of the light source acts as a spatial or temporal gate that 
allows the reference arm to be used for ranging and slice-selection. 
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Signals that are not within the temporal gate are detected as a mean 
background intensity that can be filtered out by an AC-coupled detector. 

The source 16 may include super luminescent diodes centered 
near Xo = 1.3jim so as to maximize light penetration in the tissue sample. 
5 However, any broadband light source can be used. The light from the 

source is input to a single mode optical fiber wherein a fiber optic 50:50 
beam splitter 18 directs half of the light to the reference arm 12 and half 
of the light to the sample arm 10. In both arms 10 and 12, light exiting 
the optical fibers is collimated by Faraday rotators. The collimated light 

10 in the reference arm is reflected back into the optical fiber by a retro- 

reflective prism 21 that is mounted either on a galvanometer or on a 
micro-translation stage. Translation of the retro-reflector 21 at a known 
velocity provides axial scanning as well as a means to modulate the 
OCT signal at a known, Doppler-induced carrier frequency f re f. Light in 

15 the sample arm 10 is steered by a galvanometer-based xy-scanner and is 

focused onto the tissue sample by fixed achromats. 

Light reflected by the reference reflector 21 and light reflected 
by the tissue sample 14 travels back to the beam splitter 18 and from 
there to a high sensitivity photoreceiver or detector 20. The 

20 photoreceiver or detector 20 converts the received optical signals to 

electrical signals. Within the detector 20 the signal is amplified and 
filtered. For demodulation based OCT, the signal from the detector 20 is 
processed by a lock-in amplifier tuned to the Doppler-induced carrier 
frequency of the signal. A time constant of 30^is was found to be 

25 optimal for OCT signals. Amplitude and phase outputs are then sent to a 

digital acquisition board. For A-scan based OCT processing, the signals 
from the detector 20 are pre-amplified, filtered by an analog 8-pole low 
pass Bessel filter and post amplified. The amplified A-scan signal is 
then sent to the digital acquisition board. Data acquisition may be 

30 performed with a 16 byte PCI board from Gage CompuScope 1610. For 
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demodulation based OCT, amplitude and phase from the lock-in 
amplifier are acquired through separate channels. For the A-scan based 
processing, the amplified signal from the detector 20 is acquired through 
a single channel. In both of these modes, the position feedback signal 
5 from the reference arm is passed to the trigger input of the data 

acquisition board. 

The signal from the detector 20 is digitized by an analog to 
digital converter of the data acquisition board. The digitized signal is 
then processed with a bi-directional filter. This type of filtering offers 
10 zero phase lag by first passing the signal forward through the filter and 

then backward therethrough, canceling out phase delays. For each pass, 
a third order elliptical bandpass filter is used. The filtered digital A-scan 
data is then processed in the time domain with sinusoidal models as 
follows. 

15 In accordance with the present invention, it is recognized that the 

OCT signal is locally a sinusoid. That is, over a small portion of a cycle 
of the OCT signal, the signal is sinusoidal. A non-linear sinusoidal 
model for a portion of the OCT signal is: 

7(0 = 4 sin (2^ + ^) 

20 where A is the amplitude of the envelope, I is the intensity of the OCT 

signal, fo is the frequency of the signal, and </> 0 is the phase lag. For a 
normalized A-scan where A = 1, the sinusoidal model is in the form of: 

7(0 = sin (2^ + ^ 0 ). 
25 The frequency to be determined by the processing is related to the 

known reference frequency,^/, and flow-induced Doppler shifts, Afyj OW9 
as follows: 

fo = fref + 4/" flow 

where f re fis determined by the reference arm scanning velocity, V re f 
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The amplitude A, frequency f 0 and phase lag (f> 0 are determined by 

non-linear regression using the sinusoidal model noted above. Figure 2 
illustrates twelve data points that span one-third of the signal's period. 
In this example, the signal to noise ratio is 30dB. These digital data 
5 points are processed using non-linear regression with the sinusoidal 

model above so as to fit a sinusoid to the data. The non-linear 
regression is implemented in software as follows: 

NonlinearRegress[data, 
10 aSin[2Pift + \[Phi]], {t}, 

{{a, aO},[f,fD},{\[Phi],\[Phi]0}} 5 

RegressionReport -> {BestFitParameters, ParameterCITable}] 

The regression yields model parameters that are best fit to the twelve 
15 data points as follows with corresponding error estimates and a 95% 

confidence value. 

5.03214 sin(0.505818 + 1884980 . 
Thus, the non-linear regression determines the frequency f 0 while 
avoiding the constraints of Fourier analysis. Fig. 3 illustrates the fit of 
the sinusoid to the digital data depicted in Fig. 2. 

In another embodiment of the present invention, the sinusoidal 
model is expanded to include changes in amplitude and shear rate as 
well. The expanded sinusoidal model is as follows: 

I{t) = (A + at) sin (2/r(f 0 + ai)t + <j> 0 ) 

where the additional parameters a and a model changes in amplitude 
and the rate of change of the frequency respectively. 

The parameters or coefficients of the sinusoidal models 
discussed above are determined by non-linear regression or optimization 
with trust-region algorithms optimized for a known frequency range, by 
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seeding the algorithm with the known frequency f re f. From the 
determined parameters or coefficients of the sinusoidal model, an image 
can be obtained depicting turbulent flow in small vessels. The 
determined frequency of the signal provides information relating to the 
5 velocity of blood flow at a certain point within the tissue and the 

parameter a provides information as to the shear rate at that point for an 
OCT signal. The image can be improved by applying a median filter to 
each axial profile computed with the extended model. A median 
velocity image can be computed from only 16 B-scans. Averaging 

10 multiple B-scans increases the signal to noise ratio of the velocity image. 

Although as many as 16 B-scans can be used, small vessels are visible 
when as few as 3 B-scans are averaged. 

Fig. 4 depicts velocity as a function of depth obtained from an 
OCT signal processed in accordance with the time domain method of the 

15 present invention providing an unfiltered A-scan. A-scans obtained with 

the method of the present invention are filtered using median filtering 
with a window spanning a coherence length l c to smooth the velocity 
profile. The median filtered velocity profile is depicted in Fig. 5. By 
removing those measurements that do not converge correctly, i.e., that 

20 are equal to the bounds specified in the trust-region optimization, noise 

is further reduced and results in a velocity profile depicted in Fig. 6 
which is similar to that obtained by frequency-domain processing but 
which has a higher resolution. Processing using a normalized sinusoidal 
model yields a similar velocity profile as depicted in Fig. 7. 

25 The method of the present invention provides high quality OCT 

images while being simpler and less expensive to implement. The signal 
processing method can be used to process other types of signals that 
locally are a sinusoid where the frequency is known within a given 
range. Many modifications and variations of the present invention are 

30 possible in light of the above teachings. Thus, it is to be understood 
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that, within the scope of the appended claims, the invention may be 
practiced otherwise than as described hereinabove. 

What is claimed and desired to be secured by Letters Patent is: 



